How far is it to a sudden future singularity of pressure? 
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We discuss the constraints coming from current observations of type la supernovae on cosmo- 
logical models which allow sudden future singularities of pressure (with the scale factor and the 
energy density regular). We show that such a sudden singularity may happen in the very near 
future (e.g. within ten million years) and its prediction at the present moment of cosmic evolution 
cannot be distinguished, with current observational data, from the prediction given by the standard 
quintessence scenario of future evolution. Fortunately, sudden future singularities are characterized 
by a momentary peak of infinite tidal forces only; there is no geodesic incompletness which means 
that the evolution of the universe may eventually be continued throughout until another "more 
serious" singularity such as Big-Crunch or Big-Rip. 



PACS numbers: 98.80.Cq; 98.80.Hw 



Over the past decade observations of high-redshift 
Type la supernovae (SNIa) have provided strong evi- 
dence that the expansion of the universe is accelerating, 
driven in the standard paradigm by some form of dark 
energyP, Q. Current dataQ continue to leave open the 
possibility that dark energy exists in the form of phan- 
tom energy, which may violate all energy conditions 
the null (gc 2 + p > 0), weak (gc 2 > and gc 2 + p > 0), 
strong (gc 2 + p > and gc 2 + 3p > 0), and dominant 
energy (gc 2 > 0, —gc 2 < p < gc 2 ) conditions (where 
c is the speed of light, g is the mass density in kgm~ 3 
and p is the pressure). Phantom matter may dominate 
the universe in the future and drive it towards a Big-Rip 
(BR) singularity in which all matter will be dissociated 
by gravity [J]. This is dramatically different from the 
standard picture of future cosmic evolution which sug- 
gests an asymptotically empty de-Sitter state driven by 
the cosmological constant or quintessence [f| and leading 
to the violation of the strong energy condition only. 

Phantom-driven scenarios have encouraged the study 
of other exotic possibilities for the future evolution of 
the universe. One of these possibilities appears in those 
models which do not assume any explicit form for the 
equation of state p = p(g), leaving the evolution of the 
energy density and pressure unconstrained. This free- 
dom may result in a so-called sudden future singularity 
(SFS) of pressure Q which violates only the dominant 
energy condition. The nature of a sudden future singu- 
larity is different from that of a standard Big-Bang (BB) 
singularity, and also from a Big-Rip singularity, in that it 
does not exhibit geodesic incompletness and the cosmic 
evolution may eventually be extended beyond it @, S]- 
The only physical characteristic of these singularities is 



a momentarily infinite peak of the tidal forces in the uni- 
verse. In more general models this peak may also appear 
in the derivatives of the tidal forces. It is interesting to 
note that these types of singularity are in a way simi- 
lar to yet another type, which were termed finite density 
singularities However, the crucial difference is that 
finite density singularities occur as singularities in space 
rather than in time, which means that even at the present 
moment of cosmic evolution they could exist somewhere 
in the Universe [13]. We will not discuss in detail finite 
density singularities in this paper since they basically ap- 
pear in cosmological models without homogeneity. On 
the other hand, it is worth mentioning that the sudden 
future singularities are quite generic since they may arise 
in both homogeneous and inhomogeneous hj mod- 
els of the universe. 

In order to obtain a sudden future singularity consider 
the simple framework of an Einstein- Friedmann cosmol- 
ogy governed by the standard field equations 
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where the energy-momentum conservation law 

e = -3-fe+4) . ( 3 ) 



is trivially fulfilled due to the Bianchi identity. Here a(t) 
is the scale factor, G is the gravitational constant, and 
the curvature index k = 0, ±1. What is crucial to ob- 
tain a sudden future singularity is that no link between 
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the energy density and pressure (the equation of state) is 
specified. This allows us to integrate ([3]) only by quadra- 
tures as 
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Of course ((IJ reduces to the standard expression for 
energy conservation, ga^ w+1 ^ — const., provided a 
barotropic equation of state, p = wgc 2 for constant w, 
is assumed. (The condition for phantom models, for ex- 
ample, is w < —1). 

From equations HJ)-© one can easily see that a pres- 
sure singularity p — ► =Foo occurs when the acceleration 
a — > ±oo, no matter that the value of the energy density 
g and the scale factor a(t) are regular. Since in that case 
| p |> g, it is clear that the dominant energy condition 
is violated. This condition can be achieved if the scale 
factor takes the form [6| 

a(t)=a s [l + (l-S)y m -5(l-y) n ] , y = f (5) 

with the appropriate choice of the constants 8, t s ,a s , to, n. 
Moreover, we can see that the r-th derivative of the scale 
factor ([5]) is given by 



W 



to(to — l)...(m — r + 1) 



+ ( _ ir ig"("-l)-("-r + l) {1 _ y y 



(0) 



and is related to the appropriate pressure derivative 
p(r-2) xhus, in general, it is possible that one has a 
pressure derivative p^ r ~ 2 ^ singularity which accompanies 
the blow-up of the r-th derivative of the scale factor cS r \ 
Observationally this could be manifested in, for example, 
the blow-up of the characteristics known as statefinders , 
such as jerk, snap etc. [3]. The pressure derivative sin- 
gularity p( r ~ 2 } appears when 



1 < n < r 
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and for any r > 3 it fulfills all energy conditions. These 
singularities are called generalized sudden future singu- 
larities (GSFS) and are possible, for example, in theories 
with higher-order curvature quantum corrections [l3j ]. 

Let us now return to the case of r = 2, for which 1 < 
n < 2 and we obtain sudden future singularity models 
of pressure (and obviously all of its higher derivatives) 
which lead to violation of the dominant energy condition. 
In such models, expressed in terms of the scale factor (0, 
the evolution begins with the standard BB singularity at 
t = for a = 0, and finishes at SFS for t = t s where a = 
a s = a(t s ) is a constant. (Note that we have changed the 



original parametrization of Ref. Q for the scale factor 
© using A = 5a s ). 

The standard Friedmann limit (i.e. models without an 
SFS) of (O is achieved when 8^0; hence 8 becomes 
the "non-standardicity" parameter of SFS models. Ad- 
ditionally, notwithstanding Ref. [j| and in agreement 
with the field equations H])-©, we assume that 8 can be 
both positive and negative leading to a deceleration or 
an acceleration (cf. ©) of the universe, respectively. 

It is important to our discussion that the asymptotic 
behaviour of the scale factor © close to the BB singu- 
larity at t — is given by a simple power-law <zbb = y m , 
simulating the behaviour of flat k — barotropic fluid 
models with to = 2/[3(w + l)] . This allows us to preserve 
all the standard observed characteristics of early universe 
cosmology - such as the cosmic microwave background, 
density perturbations, nucleosynthesis etc. - provided we 
choose an appropriate value of to. On the other hand, 
close to an SFS the asymptotic behaviour of the scale 
factor is non-standard, asFS = o-s [1 — S (1 — y) n ], show- 
ing that asFS = a s for t — t s (i.e. y = 1) at the SFS. 
Notice that one does not violate the energy conditions if 
the parameter to lies in the range 



0<to<1 (w>-l/S), 



(8) 



This range of values is, in fact, equivalent to a standard 
(neither quintessence-like nor phantom-like) evolution of 
the universe. However, with no adverse impact on the 
field equations H])-©, one could also extend the val- 
ues of to to lie in the complementary ranges to > 1 
(i.e. — 1 < w < —1/3) for quintessence, and m < (i.e. 
w < —1) for phantom, although these ranges may lead 
to violation of the strong and weak energy conditions 
respectively. 

We will next calculate the luminosity distance as a 
function of redshift, and hence the redshift-magnitude 
relation, for SFS models. This will allow us to estab- 
lish whether these models are a realistic possibility for 
the future evolution of the universe, and more specifi- 
cally whether current cosmological observations of high 
redshift supernovae are consistent with values of the con- 
stant n in the range 1 < n < 2, as required in order that 
the scale factor will display an SFS (or, more generally, 
a GSFS for r — 1 < n < r). We will then explore the 
range of values for the other SFS model parameters which 
are consistent with current observational constraints on 
standard cosmology, and thus determine limits on how 
far into the future an SFS might occur. In fact, as we 
will see below, we need to consider only two further pa- 
rameters: 8 and yo = to/ts, where to is the current age 
of the Universe in the SFS model. Notice that, in view of 
(j8]), it is reasonable to take to = 2/3 as for the standard 
dust-dominated evolution. This implies that, at early 
times, our SFS model reduces to the Einstein-de-Sitter 



universe. 
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We proceed within the framework of Friedmann cos- 
mology, and consider an observer located at r = at 
coordinate time t = t . The observer receives a light ray 
emitted at r = r± at coordinate time t = t\. We then 
have a standard null geodesic equation 



dr r° cdt 
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with the scale factor a(t) given by ([5|). Using ([5]) again, 
the redshift is given by 

_ a(t ) _ 6+(l-5)yF-6(l-y ) n 

where yo = y(to) and y\ = y{t\). The luminosity distance 
is defined as 
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Neglecting extinction and k— corrections, the observed 
and absolute magnitudes of a source at redshift z and 
luminosity distance Dl are related by 

m(«)=M-51og 10 ff + 25 + 51og 1 oX>L(«), (12) 

which, with the help of the equation ([5]), (|10p and (jTTj) . 
allows a redshift-magnitude relation for SFS cosmologi- 
cal models to be constructed. It is obvious that equation 
([9]) has to be integrated numerically in order to estab- 
lish the relation between to and t\, which can then be 
inserted into (fTP)) and (fTTj) to constrain the SFS model 
parameters. As a first step we determine the dependence 
on the SFS model parameters of the Hubble law, which 
replaces equation (| 1 2[) for z m 0, i.e. cz as HqDl, where 
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is the present value of the Hubble parameter, which we 
can take as 72kms~ 1 Mpc~ 1 [l|. 

Similarly we could derive an expression, in terms of 
the SFS model parameters, for the deceleration param- 
eter qo = — (da/a 2 )o- However, in order to search the 
parameter space for models which are admissible by cur- 
rent observations, we write the product of H and q as 



qoH 
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t m(m - 1)(1 - 6)y r o n ~ 2 - 6n(n - 1) (1 - yo)' 
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In order to obtain an accelerated universe at the present 
moment of the evolution, this product should be nega- 
tive. Fig. [T] shows an example plot of the product Hoqo 



as a function of S and yo, with the other parameters fixed 
at 77i — 2/3, n = 1.9993, t = 13.2457 Gyr. From the plot 
we see that there are large regions of the parameter space 
which admit cosmic acceleration. We have explored the 
parameter space further with various configurations of 
in, n, S, yo, to, qo, and Ho, and obtained the general con- 
clusion that there is a large class of SFS models which 
are compatible with current acceleration. 
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FIG. 1: Parameter space (Hoqo,8,ya) for fixed values of 
77i = 2/3, n = 1.9993, to = 13.3547 Gyr of the sudden future 
singularity models. There are large regions of the parameter 
space which admit cosmic acceleration. 



Out of these admissible models we then searched for 
those which are compatible with the redshift-magnitude 
relation (JT3J) observed for recent SNIa data Q , and hence 
with the derived parameters of the standard 'Concor- 
dance cosmology' (CC). We were able to identify SFS 
models that are in remarkably tight agreement with cur- 
rent SNIa data. As an illustrative example Fig. [2] shows 
luminosity distance as a function of redshift for the CC 
model with Ho — 72kms -1 Mpc~ 1 , fl m o = 0.26 and 
^ao = 0.74, and an SFS model with parameters m = 2/3, 
yo = 0.99936, S = -0.471, n = 1.9999. We see that the 
SFS model mimics the CC model very closely over a wide 
range of redshifts. In particular, it is clear that recent 
SNIa data from the Tonry at al. 'Gold' sample [l[ and 
SNLS sample cannot yet discriminate between the CC 
and SFS models. 

Taking the current age of the universe in the SFS model 
to be equal to the age of the CC model, i.e. to = 13.6Gyr, 
we find that the time to the sudden singularity is t s — to ~ 
8.7Myr, which is amazingly close to the present epoch. 
In that context there is no wonder that these singularities 
are called "sudden" . We have also checked that the larger 
the value of r in ([7]) the later in future a GSFS appears. 
It means that the strongest of these singularities which 
violates the dominant energy condition (i.e. an SFS) is 
more likely to become reality. 

Our remark about the effect of the sudden pressure 
singularity seems in agreement with the result of Ref. 
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log, (redshift) 

FIG. 2: The distance modulus fiL = m — M for the con- 
cordance cosmology (CC) model with Ho = 72kms _1 Mpc _1 , 
fimo = 0.26, Qao = 0.74 (dashed curve) and sudden fu- 
ture singularity (SFS) model for m = 2/3, n = 1.9999,5 = 
-0.471, y = 0.99936 (solid curve). Also shown are the 'Gold' 
(open circles) and SNLS (filled circles) SNIa data. Taking the 
age of the SFS model to be equal to that of the CC model, 
i.e. to = 13.6 Gyr, one finds that an SFS is possible in only 
8.7 million years. 



[16| which showed that the dominant energy condition is 
now violated and that it became violated quite recently 
(at redshift z ~ 0.2). Of course this violation may also 
be due to phantom energy 

In conclusion, we have shown that a sudden future 
singularity may happen in the comparatively near fu- 
ture (e.g. within ten million years) and its prediction 
at the present moment of cosmic evolution cannot be 
distinguished, with current observational data, from the 
prediction given by the standard quintessence scenario 
of future evolution in the Concordance Model. Fortu- 
nately, sudden future singularities are characterized by a 
momentary peak of infinite tidal forces only; there is no 
geodesic incomplctness which means that the evolution 
of the universe may eventually be continued beyond the 
SFS until another "more serious" singularity such as a 
Big-Crunch or a Big- Rip. One could then consider, more 
generally, a scale factor of the form 0, HH 

o(t) = A+[(a s -A)-D(t r -t s ) p -Et°}y m (15) 
- (A + DtP) (1 - y) n + D(t r - t s yf + Et° s y° , 

where the constants m, o,p, A, D, E are chosen so that 
the universe begins with a Big-Bang at t = where a = 0, 
next faces a sudden future singularity at t = t g where 
a(t s ) — a s , and then eventually continues to a Big- Rip 
at t = t r where a(t r ) — ► oo. All of the matter sources 
may be involved since the constants in (|15[) can be taken 
as: < to < 1 (quintessence), p < (phantom), and 
o > 1 (standard positive matter pressure). 

Whether the universe will end in a Big-Rip or a Big- 
Crunch is an open question. Moreover, unlike a sudden 



future singularity, both a Big-Rip and Big-Crunch sin- 
gularity would represent the real end of the universe. 
Fortunately, as was shown in Refs. 0, [13], a Big- Rip 
singularity is not possible in the very near future: in or- 
der to reach it one must wait about the same time as 
the current age of the universe. Apart from that, it is 
still possible to avoid it due to a negative tension brane 
contribution in a turnaround cyclic cosmology [181 ] . 
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